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ABSTRACT 


By  adapting  the  streamline  model  developed  by  Caudle  and  LeBlanc 
for  flow  through  porous  media,  a  streamline  model  for  in-situ  combustion 
was  developed. 

Initial  tests  were  made  to  compare  the  streamline  model  for  in-situ 
combustion  with  actual  field  trials  conducted  by  oil  companies.    The  results 
were  extremely  close.    As  a  result,   in-situ  combustion  experiments  may 
be  examined  and  optimized  to  obtain  the  best  possible  results  prior  to  im- 
plementation in  the  field. 

Conducting  numerous  tests  by  changing  assumed  values  indicated 
that  areas  in  the  process  must  be  re-examined  in  order  to  determine  the 
economic  feasibility  of  using  in-situ  combustion  to  recover  the  residual 
oil  remaining  in  the  nation's  reservoirs,  kerogen  from  oil  shale,  and  the 
hydrocarbons  from  tar  sands. 
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CHAPTER    I 
INTRODUCTION 


In  1972  the  United  States  became  aware  of  an  energy  crisis.    The 
crisis  took  the  form  of  shortages  of  fuel  for  production  of  electricity  and 
heating  fuels  for  homes  and  businesses.    The  causes  are  many  and  varied, 
but  may  be  traced  to  one  reason  —  economics.    Since  the  start  of  oil  pro- 
duction, the  cost  per  barrel  has  been  kept  so  low  by  artificial  barriers  that 
the  return  of  the  dollar  invested  in  oil  or  gas  exploration  and  production  is 
below  the  return  of  investments  in  other  areas.    Thus  the  energy  crisis  is 
not  the  lack  of  fuel  per  se,  but  the  lack  of  economic  fuels. 

All  persons  intimately  related  with  the  oil  industry  know  that  existing 
means  of  production  may  optimistically  recover  fifty  percent  of  the  oil  in 
the  nation's  reservoirs.    Thus  one  could  conservatively  estimate  that  over 
90  billion  barrels  of  oil  still  remain  in  established  reservoirs,     since  the 
total  production  of  the  United  States  from  1920  to  1969  was  just  over  the  90 
billion  barrel  mark  (1).    In  addition,  potentially   1,781   billion  barrels  of 
kerogen  from  oil  shale  and  399  billion  barrels  of  oil  from  tar  sands  exist  in 
the  United  States  (2).    Based  on  an  estimated  rate  of  consumption  of   46  mil- 
lion barrels  per  day  by  1980,  the  above  resources  could  supply  the  United 
States  for  over  135  years. 

The  fallacy  with  the  above  conclusion  is  that  the  price  of  a  barrel 
of  oil  would  have  to  rise  drastically  to  make  production  of  the  residual  oil, 
oil  shale,   and  tar  sands  profitable  if  present  production  concepts  are  main- 
tained.   For  example,  present  methods  of  oil  production  from  oil  shale  is 
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by  mining  the  oil  shale  and  retorting  the  shale  above  ground.  Not  only  is 
the  process  environmentally  distasteful,  but  the  resulting  cost  per  barrel 
of  oil  is  between  four  and  five  dollars  (2). 

During  the  past  ten  to  fifteen  years,   one  process  has  received  exten- 
sive investigation   —   in-situ  combustion,   often  called  fire  flooding  or  under- 
ground retorting.    Through  in-situ  combustion,  residual  oil  in  existing  oil 
fields  could  be  recovered.    Synthetic  oil  from  oil  shale  and  tar  sands  could 
be  recovered  without  strip  mining,  and  possibly  at  lower  costs  than  the 
present. 

To  the  author's  knowledge  no  attempt  has  been  made  to  simulate  in- 
situ  combustion  by  use  of  the  Caudle-LeBlanc  Streamline  Model  (3).    Ex- 
pressed simply,  the  Caudle-LeBlanc  model  is  a  mathematical  representation 
of  the  movement  of  one  typical  particle  of  oil.    The  following  study  describes 
the  basis  of  the  model  and  the  modification  of  the  model  to  include  the  specific 
characteristics  of  in-situ  combustion,  with  the  concluding  sections  comparing 
the  model  results  with  the  results  of  actual  field  tests  conducted  by  various 
organizations. 


CHAPTER  II 
THE  STREAMLINE  MODEL 


The  theory  and  development  of  the  streamline  model  has  been  re- 
viewed by  many  different  authors  (3,  4).     For  this  reason  only  a  brief  review 
of  the  model  will  be  shown.     Following  this  brief  introduction  the  model 
modifications  will  be  shown. 

The  Basic  Streamline  Model 

The  basic  equations  were  developed  from  Darcy's  Law  and  included 
the  following  assumptions: 

1.  negligible  gravitational  effects  inside  the  reservoir 

2.  incompressible  homogeneous  fluids 

3.  horizontal  system,   infinite  in  extent 

4.  homogeneous,   isotropic  media 

The  most  general  form  of  Darcy's  Law  is  as  follows: 


kA     dp  ,       . 

q="  —  -w  (2-« 


It  is  also  known  that  the  fluid  flux  is  equal  to  the  flow  rate  divided  by  the 
area  through  which  it  flows,  thus  Equation  2.  1  may  be  written  as 


u     =     -i-    =     -JL  iP  (2.2) 

A  n      d2 
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In  order  to  determine  the  streamline  path  it  is  necessary  to  know  its  velocity 
in  the  x,  y,  and  z  directions.    Since  it  is  known  that  the  velocity  is  the  fluid 
flux  divided  by  the  porosity  of  the  material  through  which  the  fluid  is  flowing, 
Equation  2.  2  may  be  written  as: 


Velocity     =     v     = 


u  k        dp  (2.3) 


0  </>/x      d2 


The  second  important  equation  is  the  continuity  equation.  It  is  a  partial 
differential  equation  derived  from  the  material  balance  around  a  point  in 
three-dimensional  space.    This  equation  is  written  as  follows: 


a(pux)        a(puv)        a(puz)  9o 

+        :  +  : =     "^  (2.4) 


3x  9y  3z  dt 

When  the  volumetric  flux  components  of  Darcy's  Law,  Equation  2.2    in  the 
x,  y,  z    directions,  are  substituted  in  Equation  2.4  the  following  equation 
results: 

^  ,p^      ./£»\         _*    (,5) 

ox  \  H    ox  ]         3y  \  \i    3y    /       3z  \   \x    dz  \  3t 

Applying  the  assumptions  made  earlier,  Equation  2.5  becomes: 


92P    +  i!g.    =    0  (2.6) 


3?         3y2 


which  has  become  known  as  LaPlace's  Equation.  By  transforming  to  cylin- 
drical coordinates  and  allowing  radial  flow  around  a  point,  the  equation  may 
be  solved  for  the  flow  potential  in  a  horizontal,   isotropic,   homogeneous 


5 
porous  medium.    Restoring  the  resulting  equation  to  rectangular  coordinates 
and  allowing  for  a  multiwell  system,   Equation  2.6  becomes: 


n 
4?rkh 


P(x>y)  Pm  ■    -TTT-    7qt*»   (x-xj2  +  (y-y-)2]  (2.7) 

i  =  l         L  J 


where   p       is  the  mean  reservoir  pressure. 

By  differentiating  Equation  2.7  with  respect  to    x   and  substituting  the 
resulting  derivative  into  Equation  2.3    the  velocity  at  any  point  in  a  rectan- 
gular system  becomes 

n  /         \ 

..         v  <1  (x_xi) 

vx<x>y>    =    ~ohr  I  — -o — ■ 7  <2-8> 

2^h    1  =  1  (x-x/  +  (y-yj)2 
In  a  similar  manner,  the  velocity  in  the   y   direction  becomes 


,    >  i     v        q(y-yi) 

vy<x>y>    =    "=■   L 5 5-  (2-9> 


2^h    i=l(x-Xi)2  +  (y-yi)2 


In  order  to  follow  the  particle  along  a  particular  streamline,  the  dis- 
tance of  travel  must  be  known.    Given  a  finite  amount  of  time,  the  distance 
a  particle  travels  at  a  constant  velocity  may  readily  be  calculated  as 

ds     =     vdt        or        As     =     vAt  (2. 10) 

Assuming  a  particle  was  initially  located  at  point    (xj,  yj) ,    the  location  of 
the  same  particle  after  time    At    could  be  calculated  as 
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xi+l     =     xi   +  vxAAt 

(2.11) 

Fi+1  =   n  +  vyiAt 

By  starting  on  the  wellbore  of  the  injection  well  and  systematically 
reiterating  Equations  2.11,  a  trace  of  any  particle's  movement  may  be  made. 
Figure   1    shows  the  paths  of  particles,   from  a  single  five -spot  pattern  with  a 
central  injection  well.    The  time  to  reach  a  specific  point  may  be  determined 
by  the  summation  of  the    At's    required  to  cause  sufficient  movement  of  the 
particle  to  reach  that  specific  point. 

In  a  given  two-dimensional  fluid  system,  an  infinite  number  of  par- 
ticle paths  (streamlines)  could  be  used  to  describe  the  flow  from  a  source  to 
a  sink.    However,  to  simplify  computation,  a  representative  number  of 
streamlines  may  be  chosen  to  represent  all  fluids  emanating  from  the  source. 
A  few  of  the  streamlines  from  a  source  are  shown  on  Figure  2    (5).     By 
dividing  the  distance  between  two  adjacent  streamlines  into  two  equal  parts, 
the  streamline  may  be  said  to  represent  the  flow  of  all  particles  within  the 
dividing  lines.    A  streamtube  has  thus  been  defined  and  is  shown  on  Figure  3. 


O  Production  Well 
•  Injection  Well 


FIGURE    1     Unbounded  Five-Spot  Pattern  Shwoing  Generated 
Streamlines  from  Wells  Having  Equal  Rates 
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FIGURE    2     Typical  Streamlines  in  a  Two-Well  System 
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FIGURE    3     Typical  Flow  Tube 


CHAPTER    m 
IN -SITU  COMBUSTION  MODEL 


Prior  to  developing  the  modifications  necessary  to  convert  the  stream- 
line model  to  a  streamline  model  for  in-situ  combustion,   a  discussion  of  the 
process  of  in-situ  combustion  is  warranted. 

In-situ  combustion  is  not  a  new  process.    Recovery  by  a  combustion 
or  heat  wave  process  was  patented  in   1923.     F.  A.    Howard  was  granted  a 
patent  on  a  process  in  which  air  and  a  combustible  gas  were  pumped  into  an 
injection  well  and  ignited  (6).    The  method  involves  ignition  of  the  formation 
in  an  injection  well,   followed  by  propagation  of  a  combustion  front  through 
the  reservoir.    Combustion  is  maintained  by  the  injection  of  an  oxygen- 
containing  gas,  such  as  air,  to  react  with  reservoir  hydrocarbons.    As  the 
flame  front,  which  is  really  glowing  embers  of  the  reservoir  material,  pro- 
gresses through  the  reservoir,   oil  and  formation  water  are  vaporized,  driven 
forward  in  the  gaseous  phase  and  recondensed  in  the  cooler  part  of  the  for- 
mation. In  turn  the  condensed  fluids  push  oil  into  the  producing  wells  (7). 

Since   1923  field  experiments  have  been  conducted  by  the  Russians  in 
1935  (8)  and  others  in  Oklahoma  (9)  and  Kansas  (7),   with  the  most  recent 
test  occurring  in  the  Green  River  Oil  Shale  deposits  by  the  U.S.  Bureau  of 
Mines  (11,  12).    The  objective  of  the  early  tests  and  those  occurring  in  the 
later  years  was  to  recover  the  residual  oil  left  after  primary  drive  and 
secondary  methods  had  been  utilized.    In  the   1950's  in-situ  combustion  ex- 
periments were  conducted  on  oil  shales  with  the  intent  of  retorting  the 
shale  underground  as  an  alternative  to  mining  the  shale  and  retorting  above 
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ground.  Other  tests  were  conducted  in  highly  viscous  oils  with  the  intention 
of  reducing  the  viscosity  of  the  oil  in  order  to  make  the  oil  more  mobile  and 
thus  producible  in  a  smaller  period  of  time  (10). 

The  complex  nature  of  the  in-situ  combustion  process,   with  it's 
many  inter-related  displacement  mechanisms,   makes  exact  interpretation 
of  field  data  impossible.    This  is  understandable  considering  that  even  two- 
phase  flow  of  fluids  through  a  porous  homogeneous  medium,   which  has  been 
studied  for  many  years,   cannot  be  adequately  described.    As  a  result,  as- 
sumptions must  be  made  based  on  experience  and  the  results  of  many  trials. 
Gates  and  Ramey  (10)  conclude  that,   as  a  result  of  their  interpretation  of 
field  tests  in  California,   Figure  4  represents  the  in-situ  combustion  process, 
Allred's  (13)  beliefs  are  somewhat  different  and  describe  the  in-situ  com- 
bustion process  as  shown  on  Figure  5.    With  the  exception  of  the  gases 
shown  in  Figure  5,  both  representations  are  the  same.    This  author  has 
simplified  the  two  figures  to  enable  utilization  in  the  streamline  model  simu- 
lation.   The  in-situ  combustion  process  as  utilized  in  the  simulation  is 
shown  on  Figure  6. 

Along  with  the  composition  of  the  various  zones,  the  velocity  of  the 
fronts  is  important.    Since  the  burn  front  is  the  driving  mechanism  in  the 
process,  its  velocity  is  the  basis  for  determining  the  velocities  of  the  other 
zones.    Two  experiments  indicate  that  the  injected  air  or  gas  must  reach  a 
certain  rate  before  movement  of  the  front  will  begin.    Figure  7    represents 
the  results  of  the  aforementioned  experiments  (9,  14).     Further  experiments 
may  show  that  the  rate  of  advancement  of  the  combustion  front  depends  on 
the  oil  content  of  the  reservoir,   reservoir  porosity  or  some  other  variables. 
What  is  important  at  this  time  is  that  the  relation  between  air  flux  across 
the  burn  and  the  rate  of  movement  of  the  burn  does  exist. 
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Other  such  relations  and  assumptions  had  to  be  established  in  the 
modification  of  the  streamline  model,  but  to  ensure  their  complete  compre- 
hension they  will  be  discussed  as  they  affect  the  streamline  model  modifi- 
cation. 

In  developing  the  original  streamline  model  we  assumed  that  the  fluid 
was  incompressible.  The  effect  of  this  assumption  is  extremely  close  to  reality. 
Merrick  (18)    compared  the  results  of  flow  of  real  gases,   ideal  gases,   and 
incompressible  fluids  through  a  porous  media  and  determined  that  the  dif- 
ferences in  results  were  extremely  small.    This  assumption  will  also  be 
made  for  the  in- situ  combustion  model  of  flow  in  a  reservoir,  but  to  utilize 
the  relationship  of  the  advancement  of  the  combustion  front  and  the  injected 
air  flux,  the  injected  air  flux  must  be  in  surface  volumes.    To  properly 
change  the  injection  rate  to  surface  volumes  we  must  look  at  the  gas  laws. 
For  an  ideal  gas 

pV     =     nRT  (3.1) 

where  V  is  the  volume  of  a  given  number  of  gas  moles  n  at  absolute  pres- 
sure p  and  temperature  T.  R  is  the  universal  gas  constant.  For  a  non- 
ideal  gas  the  relation  is 

pV     =     znRT  (3.2) 

where    z    is  the  compressibility  factor  of  the  gas.    The  compressibility  is 
a  function  of  the  pressure  in  the  reservoir    pm.    Thus,  Equation  2.3    de- 
scribing an  ideal  gas  becomes 

v     =      U  Zm  (3.3) 


</>p 


m 
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where  z  is  the  compressibility  factor  of  the  fluid  injected  into  the  reser- 
voir, and  p  is  the  average  reservoir  pressure  in  atmospheres  (standard 
pressure  is  one  atmosphere)  and  the  surface  and  reservoir  temperature  are 

equal. 

Equations    2.7,   2.8,    and  2.9    must  also  be  modified  since  they 

were  developed  for  noncompressible  flux.    With  an  ideal  fluid,  the  volume 

at  standard  conditions  equals  the  volume  at  reservoir  conditions.    Since  we 

are  using  a  compressible  fluid,  the  quantity  at  surface  conditions  becomes 

pm 

qsc  zM    qres  (3.4) 

m 

Since  it  is  also  known  that  the  mobility  of  a  fluid  is  equal  to  its  permeability 
divided  by  its  viscosity,  Equations    2.7,     2.8,  and  2.9   become 


P       =       Pm    " 


zm 


m       47rhpmXs 


^qSCi  ^[(x-xi)2  +  (y-yi)2j  (3.5) 


z  1  (x"xi) 

zm 


Iqsci"; — j—, — X"  (3*6) 


sx  277h0pm  1t1      *    (x-x^  +  (y-yi) 


zm      £  (y-yj) 

vs„    =     TTul Z.qsci   ~,        71       ;        TT  (3.7) 

&y  27rh0pm  gx     i    (x-xjr  +  (y-yir 

From  the  above  equations  the  pressure  and  velocity  of  the  standard  fluid  may 
be  determined.    The  standard  fluid  in  this  thesis  will  be  the  air  injected  into 
the  reservoir.    It  is  called  the  standard  fluid  to  keep  the  process  as  general 
as  possible.    In  addition,  the  standard  fluid  has  particular  significance  in  the 
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conductivity  ratio  which  will  be  explained  in  the  coming  pages. 

To  determine  the  amount  of  production,  the  velocity  of  the  burn  is 
necessary  since,   as  will  be  shown  later,  the  movement  of  all  the  zones 
containing  oil  are  dependent  on  the  movement  of  the  burn.    Equations    3.6 
and  3.7    will  generate  the  velocity  of  the  standard  fluid,  which  in  this  case 
will  be  the  air  injected  into  the  reservoir.    Knowing  the  velocity  of  the  air, 
and  by  the  relationship  of  the  air  velocity  and  burn  front  velocity  as  expres- 
sed in  Figure  7,    the  velocity  of  the  burn  may  be  computed.     From  Figure  7 
the  velocity  of  the  burn  is 


vb       =     SLvs      +   C  (3.8) 

x  x 


vb       =     SLv.      +   C  (3.9) 

y  y 

where     SL     is  the  slope  of  the  line  and   C    is  the  line  constant.    Note  that 
at  approximately  thirty  feet  per  hour,  the  relation  between  the  burn  and  the 
standard  fluid  changes.    At  this  point  the  values  of   SL    and    C    will  be 
changed.    In  the  streamline  computer  program  these  values  are  used  as  input 
data  to  provide  as  much  flexibility  as  possible.    As  previously  mentioned, 
insufficient  testing  has  been  completed  to  prepare  exact  figures  for  the  burn 
front  movement  with  regard  to  the  rate  of  injection  of  air,  but  the  fact  that 
the  relation  does  exist  must  be  considered. 

The  initial  streamline  model  developed  in  Chapter  II  was  for  a  single 
fluid.    For  a  multifluid  condition,   such  as  a  waterflood  or  in-situ  combustion 
where  different  fluids  exist,  the  concept  of  the  conductivity  ratio  must  be 
applied.    Caudle  (15)  presents  a  detailed  explanation  of  the  concept  and  thus 
it  will  not  be  discussed  here.    The  conductivity  ratio  is  determined  as  fol- 
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lows: 

y    =  -r^  (3.10) 


Pi" 

•pP 

X* 

L 

L 

<Pk- 

•  l 

Pk) 

^s 

K 

k=l 


where    y   is  the  conductivity  ratio,    pi    is  the  pressure  at  the  source  or  in- 
jection well,    p      is  the  pressure  at  the  sink  or  production  well,    \s    is  the 
mobility  of  the  standard  fluid,  and   Xk    is  the  mobility  of  the  fluid  in  the    k 
zone  which  has  a  pressure    pk    at  its  interface  with  the  zone  preceding  the 
k    zone.      L    is  the  number  of  different  zones. 

Substituting  Equation  3.  3  in  Equations  3.  8  and  3.9  and  multiplying  by 
the  conductivity  ratio 


/  P 

vb       =    SL(vs    0—    r]+   C  (3.11) 

x  V     x      zm       / 


vbv     =    SL(vs    0—  y)   +   C  (3.12) 

y  v  Dy     zm     J 

As  a  result  the  velocities    (x  and  y    directions)  may  be  computed  for 
the  burn  front.    The  velocities  of  the  other  three  fronts  will  be  similar,  but 
will  differ  by  some  factor  due  to  the  accumulation  of  oil  in  front  of  the  burn. 
By  referring  to  Figure  8  it  may  be  seen  that  the  velocity  across  point    B 
will  yield  a  quantity   Q,     of  a  known  amount  if  the  cross-sectional  area  of 
point    B    is  known.    By  the  same  rationale  a  fluid  moving  across  point    2 
will  yield  a  quantity  of  fluid   Qq    of  known  amount  if  the  area  and  velocity 
of  point    2    are  known.    Additionally,   if  a  constant  rate  of  injection  is  estab- 
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lished,  the  quantity  of  fluid  passing  point  B  will,  in  some  way,  be  related 
to  the  quantity  of  fluid  passing  point  2.    Mathematically, 

Ob     =     vbAb  (3.13) 

Q2     =     v2A2  (3.14) 

Q2     =     F2Qb  (3.15) 
Thus,  by  substituting  Equation  3.  13  in  Equation  3.  15 

«2    =    F2vbAb  <3-16> 

Considering  the  constant  rate  of  injection,  the  following  conclusions  may  be 
made: 


Qs     =     vSbAb  (3.17) 


Qs     =     vS2A2  (3.18) 


Equating  Equations  3. 17  and  3.  18  and  solving  for    A^ , 


vsbAb 
A2     =     (3.19) 

Vs2 


Substituting  Equation  3.  19  in  Equation  3.  16, 


vsbAb 
v2     =    (^ 1        F2vbAb 
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Solving  for  the  velocity  at  point    2, 


v2     =     F2vb—  (3-21) 


By  similar  reasoning  point  3  on  Figure  8  would  have  a  velocity  of 

Vs3 

v3   =   Vbir1  <3-22) 

6  vb3 

As  shown  on  Figure  6,  the  fourth  zone  is  composed  entirely  of  gas,  and  by 
a  piston-like  displacement  the  zone  is  pushed  by  the  air  injected  into  the 
reservoir.    Thus  the  velocity  of  zone  4  is  that  of  the  standard  fluid  at  zone 
4  as  computed  by  Equations  3.6  and  3.7. 

Considering  all  the    x   and   y    components  of  the  velocities  of  zones 
2,  3   the  movement  of  these  two  points  may  be  computed  as  follows: 


v 


F2vbxvs2x 

-i?  (3.23) 


x  vc 

^x 


F2vb   vs 
I   Dy    s2y 

v9       =     J—  (3.24) 


*by 


F3VbxVs3x 


v,       =     —  (3.25) 

v  V^. 

5bx 


x  v 
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FIGURE    8     Process  Zones  in  a  Non-Linear  System 


FIGURE    9     Process  Zones  in  a  Linear  System 


F3vb   % 
y      3y 


^ 
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(3.26) 


To  compute  the  values  of    F«    and    F„  ,     a  movement  of  the  burn  must 
be  considered.    Referring  to  Figure  9,  a  movement  of  the  burn  from    B    to 
B'    will  cause  the  displacement  of  all  oil  existing  in  the  area  between    B    and 
BT   that  is  not  consumed  as  fuel  in  the  combustion  process.    Thus 


Since 


Oil  moved     =     V 


B 


0T-0 


JM 


(i.o-sgi-swi) 

(1.0 -Sg.) 


(3.27) 


dVB         (0T-0g)(l.O-fb)(l.O-Sg.-SWi) 


dt 


> 


(1.0 -Sg.) 


(3.28) 


Vi    =     A 


i^eff 


(3.29) 


Vi  dt  A  0 


eff 


> 


dVB        I  (0T-0g)(l.O-fb)(l.O-SgrSw.) 


(i-sgi) 


(3.30) 


Since 


dVB/Ai 


dt 


(3.31) 
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Equation    3.30    may  be  written  as 


(0T-0g)(l.O-fb)(l.O-Sg.-Sw.)- 
0eff(l.O-Sgi) 


(3.32) 


where    0   ~    is  the  effective  porosity  through  which  the  fluid  must  move  and 
is  computed  as 


0eff    =     0gU.O-Sg-Sw) 


(3.33) 


As  was  discussed  earlier,  the  velocity  of  the  burn  is  proportional  to  the 
velocity  of  the  air  flux  across  the  burn  front.    Thus  if  we  consider  point  2 
as  the  i-th  point, 


v 


% 


(3.34) 


Vo     v 
s2 


s«  vb 


'*2 


v, 


(3.35) 


and  Equation  3.  32  becomes 


(0T-0)(l.O-fb)(l-S      -s     ) 


W1-^' 


V         V, 

s2  vb 


(3.36) 


Comparing  Equation  3.36  with  3.21,    we  find  that    F2    is 
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(0T-0g)(l.O-fb)(l.O-Sg2-SW2) 


V1^ 


(3.37) 


By  mutual  displacement,  the  movement  of  point  2  to    2T    has  moved 
point  3  to  3'    by  an  amount  equal  to    F2-    Thus  the  movement  of  point    3    to 
3'    is 


v3'     =     F2 


vs3'vb 


sb 


(3.38) 


The  movement  of  point  3'    to    3"    caused  by  the  movement  of  the  burn  and 
not  the  movement  of  point    2  to  2f    is  given  by  Equation  3.  32  as 


Vol       = 


(0T-0g)(i.o-fb)(i.o-sg3f-sW3t)- 


w^-v 


vs3,vb 


(3.39) 


Combining  Equation  3.  38  and  3.  39  to  obtain  the  total  movement  of  point    3 
to    3"  , 


r(0T-0g)(l.Ofb)(l.O-Sg3-SW2) 


+  F, 


Vs3Vb 


% 


(3.40) 


Comparing  Equation  3.  40  to  Equation  3.  22,     F3    is  found  to  be 
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(0T-0g)(l-O-fb)(l-S      -S      ) 

F3     =     +    F2  (3.41) 

0eff(l-Sg3) 


Substituting  the  equations  developed  in  this  chapter  into  the  streamline 
model  developed  by  Caudle  and  LeBlanc  (3)  will  result  in  an  in-situ  combus- 
tion streamline  model. 


CHAPTER    IV 
EXPERIMENTAL  INVESTIGATION 


Sensitivity  to  Assumptions 

As  was  discussed  earlier,   little  is  known  about  the  in-situ  combustion 
process.    Present  technology  prevents  us  from  obtaining  cores  from  the 
reservoir  as  the  process  moves  through  it.    As  a  result,  numerous  assump- 
tions had  to  be  made.    The  first  assumption  is  the  relationship  between  the 
advancement  of  the  combustion  front  and  the  injected  air  flux.    The  effect  on 
production  and  time  of  breakthrough  by  changes  in  this  relationship  must  be 
determined.    The  second  assumption  is  the  amount  of  reservoir  hydrocarbon 
consumed  in  the  propagation  of  the  combustion  front.    As  may  be  seen  in 
Equations  3.36  and  4.40,  the  velocity  of  zones  2    and  3    are  directly  re- 
lated to  the  amount  of  reservoir  hydrocarbon  consumed.    The  third  variable 
is  the  mobility  of  the  various  zones  within  the  process. 

In  order  to  determine  the  sensitivity  of  the  streamline  model  to  vari- 
ations in  these  assumptions,   a  basic  test  was  conducted  with  the  following 
data. 

1.  The  rate  of  advancement  of  the  combustion  front  in  relation  to 
air  flux  is  in  accordance  with  the  curve  in  Figure  7  established  by  Martin 
(14)  because  it  shows  a  realistic  slowing  of  the  combustion  front  as  the  air 
flux  declines. 

2.  The  saturation  of  gas,  oil,  and  water  in  each  zone  is 
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Zone  1:         Sg     =    1.00 

Zone  2:         Sg  =0.33 

S0  =  0.34 

Sw  =  0.33 

Zone  3:         Sg  =  0.50 

S0  =  0.30 

Sw  =  0.20 

Zone  4:         Sg  =  1.00 

3.  The  rate  of  injection  was  established  at  600  scf/hr. 

4.  The  mobility  of  each  zone  was  assumed  to  be  25,    11,    12,   25,   25 
darcy/centipoise  for  Zones  1,  2,  3,  4,   and 5    respectively.    The  wisdom  of 
using  these  mobility  ratios  in  light  of  the  saturation  and  make-up  of  each 
zone  may  be  dubious;  however,  tests  were  conducted  to  determine  the  effect 
of  changes  in  zone  one  mobility  only.    The  others  were  not  examined  due  to 
reasons  to  be  explained  later. 

5.  The  amount  of  oil  consumed  in  the  propagation  of  the  combustion 
front  was  assumed  to  be  ten  percent  of  the  total  hydrocarbons  in  the  reser- 
voir. 

6.  The  average  amount  of  oil  in  the  reservoir  was  assumed  to  be 
27  gallons  of  oil  per  ton  of  reservoir  material. 

7.  Pattern  utilized  was  an  unbounded  inverted  five  spot,    100  feet 
square . 

The  above  data  were  incorporated  into  the  computer  program  shown 
in  Appendix  B.    The  plot  of  the  streamlines  is  shown  in  Figure    10,     and  the 
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FIGURE    10     Streamlines  in  Basic  Reservoir 
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production  data  are  shown  in  Table  I   with  the  location  of  the  oil  bank  at 
time  of  breakthrough  in  Figure  11.    In  subsequent  tests  the  location  of  the 
combustion  front  will  be  compared  with  the  basic  model  plot  in  order  to 
show  the  amount  of  the  reservoir  swept  by  the  burn. 

The  relationship  between  the  advancement  of  the  combustion  front 
and  the  air  flux  is  fundamental  to  successful  simulation  of  in-situ  combus- 
tion.   The  two  curves  showing  this  relationship  in  Figure  7    have  been  pub- 
lished as  a  result  of  experimental  studies  (9,  14).    As  a  result,  these  curves 
must  be  assumed  to  include  such  effects  as  oxidation  potential  of  the  oil, 
heat  conduction  and  heat  convection.    The  basic  program  was  computed 
using  the  relation  developed  by  Martin  (14).    To  determine  the  significance 
of  the  difference  in  the  Martin  and  Moss  (9)  curves,   the  second  test  was  to 
use  all  the  data  from  the  basic  test  with  the  Moss  curve.    Figure   12    com- 
pares the  results  of  the  tests,  where  the  basic  data  using  the  Martin  curve 
is  a  solid  line  and  the  second  test  using  the  Moss  curve  is  a  broken  line. 
From  the  start  of  ignition  to  129  hours,  the  fronts  are  propagated  radially. 
The  positions  of  the  fronts  shown  in  the  upper  right  quadrant  is  the  position 
of  the  combustion  fronts  at  breakthrough.    The  number  of  hours  shown  in 
that  quadrant  is  the  time  of  breakthrough  of  the  indicated  front.    Table  II 
shows  the  production  data  for  the  second  test.    The  difference  in  the  break- 
though  time  and  the  total  production  may  be  explained  by  review  of  Figure  7. 
At  high  air  flux  rates,  the  front  of  the  basic  test  moves  at  a  faster  rate. 
Production  data  confirms  this  in  that  the  initial  oil  production  is    97  percent 
greater  in  the  basic  test,   indicating  a  larger  area  swept  in  the  same  amount 
of  time.    As  the  air  flux  rate  drops  due  to  the  expanding  surface  area  of  the 
burn  front,  the  Martin  curve  declines  more  rapidly  than  the  Moss  curve. 
After  629  hours  the  combustion  front  of  the  basic  test  has  reached  the  pro- 


TABLE    I 
RECOVERY  PERFORMANCE  FOR  BASIC  TEST 


Gas  Produced 

Oil  Produced 

Water  Produced 

Time 

(ft3) 

(bbls) 

(bbls) 

(hrs) 

0.0 

0.0 

0.0 

4.2 

0.0 

0.0 

0.0 

45.2 

0.0 

0.0 

0.0 

87.5 

53.3 

3.3 

2.4 

129.0 

143.0 

10.0 

8.3 

171.0 

252.0 

18.2 

15.8 

212.0 

376.0 

27.8 

24.7 

254.0 

506.0 

38.2 

34.4 

296.0 

649.0 

49.4 

44.8 

337.0 

779.0 

59.6 

54.6 

379.0 

948.0 

72.6 

66.4 

421.0 

1070.0 

83.0 

76.4 

462.0 

1190.0 

92.9 

86.0 

504.0 

1360.0 

106.0 

98.0 

546.0 

1480.0 

115.0 

106.0 

587.0 

1580.0 

122.0 

114.0 

629.0 
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FIGURE    11     Location  of  Oil  Zone  in  Basic  Test  at  Break- 
through 
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Basic  Test  (Martin  Relationship) 
Second  Test  (Moss  Relationship) 


45.8  hrs 


1130.0  hrs    \\ 
629.0  hrs 


J O 


129.0  hrs 


FIGURE    12 


Combustion  Front  Location  Using  Martin 
and  Moss  Relationship  Between  Injected 
Air  Flux  and  Rate  of  Frontal  Advance  (14,  9) 


TABLE    II 
RECOVERY  DATA  FOR  SECOND  TEST 


Gas  Produced 

Oil  Produced 

Water  Produced 

Time 

(ft3) 

(bbls) 

(bbls) 

(hrs) 

0.0 

0.0 

0.0 

4.2 

0.0 

0.0 

0.0 

45.8 

0.0 

0.0 

0.0 

87.5 

56.8 

0.  1 

0.1 

129.0 

145.0 

2.8 

1.9 

171.0 

242.0 

7.8 

5.2 

212.0 

356.0 

14.2 

9.4 

254.0 

481.0 

21.8 

14.7 

296.0 

606.0 

30.5 

21.6 

337.0 

734.0 

40.4 

30.  1 

379.0 

864.0 

50.6 

39.0 

421.0 

1000.0 

61.7 

49.0 

462.0 

1140.0 

73.1 

59.4 

504.0 

1280.0 

85.3 

70.6 

546.0 

1430.0 

98.2 

82.5 

587.0 

1580.0 

112.0 

95.3 

629.0 

1740.0 

126.0 

108.0 

671.0 

1910.0 

140.0 

121.0 

712.0 

2070.0 

154.0 

134.0 

754.0 

2240.0 

169.0 

149.0 

796.0 

2400.0 

184.0 

163.0 

837.0 

2550.0 

198.0 

177.0 

879.0 

2670.0 

211.0 

189.0 

921.0 

2800.0 

224.0 

202.0 

962.0 

2920.0 

235.0 

213.0 

1000.0 
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duction  well  and  the  streamlines  midway  between  two  production  wells  have 
stopped.  The  reason  for  these  streamlines  stopping  may  be  found  in  Equa- 
tion 3. 10  which  has  been  expanded  below  to  include  all  zones  considered  in 
this  thesis. 


y    = 


Pi-P 


p 


pi"Pb         Pb"P2  P2-P3        P3-P4  P4-P1 

+    — +    — : +    — + 


A 


X, 


x. 


A. 


X, 


A, 


(4.1) 


As  the  streamlines  between  two  production  wells  move  out  from  the  injection 
well,  the  pressure  at  those  points  decreases  which  in  turn  decreases  the 
conductivity  ratio.    Since  the  conductivity  ratio  is  directly  proportional  to 
the  velocity  of  the  burn,  as  shown  in  Equation  3.11,     a  decrease  in  the  con- 
ductivity ratio  will  cause  the  combustion  front  to  slow  and  eventually  stop. 
The  same  would  be  true  for  the  streamlines  reaching  the  production  wells, 
except  that,  as  each  zone  is  produced,   its  term  in  the  denominator  of  Equa- 
tion 4.1    is  dropped  thus  causing  an  increase  in  the  conductivity  ratio,    The 
slowing  and  stopping  of  the  streamlines  between  the  production  wells  causes 
the  oil  front  ahead  of  the  burn  to  stop  before  reaching  the  production  wells. 
The  next  area  of  investigation  was  to  determine  the  effect  of  varia- 
tions in  the  amount  of  hydrocarbon  consumed  in  the  propagation  of  the  com- 
bustion front.    Equations  3.38    and    3.42    show  the  relationship  between  the 
velocity  of  the  different  fronts  to  the  amount  of  reservoir  hydrocarbon  con- 
sumed.   As  the  amount  of  reservoir  hydrocarbons  consumed  is  increased, 
the  velocity  of  each  front  will  be  decreased.    Figure   13    shows  that  as  the 
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FIGURE    13     Effect  of  Variation  in  Hydrocarbon  Consumption 
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amount  of  hydrocarbon  consumed  is  increased  from  10  percent  to  60  percent 
the  amount  of  production  and  the  time  required  for  breakthrough  increase. 
When  the  consumption  rate  is  set  at  70  percent  of  the  reservoir  hydrocarbon, 
the  total  production  and  time  of  breakthrough  decrease. 

The  results  indicate  that,  as  the  velocity  of  the  combustion  front  de- 
creases, the  amount  of  production  also  decreases.    Confirmation  of  this 
effect  should  be  obtained  when  the  injection  rates  are  varied.    However, 
prior  to  studying  the  effect  of  variation  in  injection  rates,  the  last  assump- 
tion should  be  discussed. 

The  third  assumption  made  in  the  utilization  of  the  streamline  model 
is  the  mobility  of  zone  1.    As  the  combustion  front  moves  across  the  reser- 
voir, the  water  and  oil  have  been  driven  forward  or  been  used  as  fuel  for 
the  combustion  process.    In  shale  the  clay  would  shrink  much  as  a  brick  does 
during  firing.    Logically,   the  permeability  of  the  zone  behind  the  combustion 
front  should  increase  and  thereby  enable  the  injected  gas  to  pass  through  the 
zone  more  easily.    The  mobility  of  this  zone  is 

X     =     —  (4.2) 

where   k    is  the  permeability  of  the  material  through  which  a  gas  or  liquid  of 
a  specific  viscosity    (/i)    is  passing.    With  an  increase  in  permeability,  the 
mobility  should  also  increase,   assuming  the  viscosity  of  the  gas  or  liquid 
does  not  change.    The  mobility  of  the  first  zone  has  the  greatest  effect  on  the 
velocity  of  the  combustion  front  since  it  is  the  first  zone  which  contains  the 
standard  fluid.    Equation  3.  10    shows  the  rationale  for  the  above  statement. 
Figure  14    shows  the  results  for  four  trials  using  a  mobility  of  the  first  zone 
of  25,   50,   75,    100  and  200.    As  may  be  seen,  the  mobility  of  zone  one  has 
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little  effect  on  the  time  to  breakthrough,  but  does  affect  the  amount  of  oil 
produced.     For  the  streamline  approaching  the  production  well  fastest,  the 
conductivity  ratio  changes  as  the  most  advanced  zone  is  produced.    Equation 
4.  1   was  expanded  to  include  all  zones  considered  in  this  paper.    As  may  be 
seen,   if  the  leading  zone  is  produced,  the  conductivity  ratio  expression  de- 
creases by  one  term  in  the  denominator  since  it  no  longer  exists  in  the 
reservoir  for  that  particular  streamline.    Since  the  mobility  of  zone  one  is 
the  same  as  the  standard,  the  conductivity  ratio  approaches  unity  as  each 
zone  is  produced,   regardless  of  the  value  used  as  the  mobility  of  zone  one. 
Thus  the  mobility  of  zone  1    has  little  effect  on  the  streamlines  that  travel 
the  shortest  distance  to  the  production  wells.    The  effect  on  the  production 
may  also  be  explained  by  Equation  4.1    in  the  same  way  as  the  stopping  of 
the  burn  front  was  determined.     However,  the  larger  the  mobility  of  the  first 
zone,  the  larger  the  value  of  the  denominator  and  the  smaller  the  conductivity 
ratio.    Thus,  the  greater  the  value  of  the  mobility  of  zone  1  the  shorter  the 
distance  from  the  injection  well  to  the  final  position  of  the  burn  front. 

The  mobility  of  the  other  zones  may  also  change  during  the  in-situ 
combustion  process,  but  their  effect  should  not  be  so  significant  as  a  change 
in  zone  1   mobility  since  they  are  filled  with  a  greater  percentage  of  liquid 
and  the  temperature  decreases  rapidly  between  combustion  front  and  zone  2. 

Injection  Rate  Variation 

Additional  tests  were  conducted  to  determine  the  effect  of  increased 
injection  rates  on  the  total  production  and  the  time  of  breakthrough  of  the 
combustion  front.     For  these  tests,   the  relationship  between  injected  air 
flux  and  the  advancement  of  the  combustion  front  was  as  developed  by  Martin. 
The  results  of  these  tests  are  shown  on  Figures  15  and  16.    These  results 
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FIGURE    16     Effect  on  Time  of  Breakthrough  by  Variation 
of  Injected  Air  Flux 
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were  rather  interesting  in  that  at  higher  injection  rates  the  total  production 
was  lower.    The  explanation  of  this  may  be  found  in  Figure  7.    At  high  in- 
jected air  flux  rates  the  combustion  front  moves  rapidly  between  the  source 
and  sink.    The  movement  is  much  faster  than  the  case  of  lower  injection 
rates  and  results  in  an  early  breakthrough.    As  a  result,  the  streamlines 
off  the  direct  line  between  the  injection  and  production  wells  are  not  able  to 
move  very  far  from  the  injection  well.    Thus,   a  smaller  amount  of  oil  is 
produced.    The  opposite  action  occurs  at  lower  injection  rates.    The  stream- 
lines on  the  direct  line  between  the  injection  well  and  production  well  move 
rapidly,  but  do  not  breakthrough  at  the  production  well  before  the  other 
streamlines  have  moved  a  considerable  distance  from  the  injection  well. 
The  graphic  evidence  of  this  action  is  shown  in  Figure    17    where  the  broken 
line  represents  the  location  of  the  combustion  front  for  an  injected  air  flux 
of   4000  scf/hr,     and  the  solid  line  represents  the  combustion  front  for  an 
injected  air  flux  of    600  scf/hr    at  breakthrough. 

The  results  of  the  tests  on  the  effect  of  increased  injection  rates  in- 
dicate that  economic  considerations  must  be  made  to  determine  the  maximum 
injection  rate  versus  the  time  and  amount  of  oil  recovered. 

Bounded  System 

Heretofore,   the  reservoir  has  been  considered  infinite  in  length  and 
width.    Now  the  reservoir  will  be  examined  with  a  boundary  existing  just 
outside  the  production  wells.    The  pattern  boundary  is  shown  in  Figure   18. 

In  previous  tests  we  noted  the  decrease  in  oil  production  due  to  the 
fingering  effect  and  rapid  expansion  of  the  combustion  front,  thereby  forcing 
oil  out  of  the  pattern  area  thus  making  it  non-recoverable.    As  is  shown  on 
Figure   11,     a  majority  of  the  oil  within  the  pattern  has  been  pushed  away 
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Injection  Rate  of  4,  000  ft3/hr 
Injection  Rate  of  600  ft3/hr 


FIGURE    17 


Combustion  Front  Location  at  Breakthrough 
with  Variation  in  Injected  Air  Flux 
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FIGURE    18     Typical  Pattern  Showing  Imaging  Technique 
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from  the  injection  well  and  in  fact  away  from  the  production  wells.    In  order 
to  close  the  system,  the  imaging  technique  will  be  utilized.    The  imaging 
technique  is  a  well-known  method  for  confining  streamlines  inside  a  de- 
sired boundary.    The  technique  is  described  in  detail  by  Collins  (13).    In 
cases  of  constant  reservoir  thickness  only  the  flow  rate  and  location  of  each 
well  is  necessary. 

The  pattern  used,   including  the  image  wells  is  shown  on  Figure  18. 
The  location  of  the  burn  front  at  breakthrough  is  shown  on  Figure  19.    Com- 
paring the  bounded  system  which  has  the  same  properties,   injection  and 
production  rates  as  the  basic  model  shows  that  the  start  of  production  and 
breakthrough  occur  fastest  in  the  bounded  system,  and  that  the  production 
rate  of  the  bounded  system  lags  that  of  the  unbounded  system.    This  could 
be  expected  since  the  imaging  technique  prevents  flow  across  the  established 
boundaries  which  are  shown  by  the  broken  line  in  Figure  18.    As  a  result, 
the  oil  laden  zones  are  directed  towards  the  production  wells;  however,   in 
the  unbounded  system,  the  zones  ahead  of  the  burn  take  a  more  circuitous 
route  as  shown  in  Figure  1    and  sweep  part  of  the  reservoir  outside  the  area 
of  the  test. 

Reviewing  the  production  data  in  Table  in  for  the  bounded  system 
shows  the  start  of  production  and  breakthrough  of  the  combustion  front  occur 
earlier  than  in  the  basic  unbounded  system.    The  bounded  reservoir  pro- 
duces about  one-third  as  much  oil  as  the  unbounded.    The  importance  of  the 
existence  of  natural  boundaries  is  therefore  apparent.    The  existence  of 
faults  close  to  the  pattern  could  cause  similar  reductions  in  production  as 
the  artificial  boundary  utilized  in  the  above  test.     Fortunately,   Lin  (16)    has 
developed  a  technique  for  construction  of  boundaries  for  the  streamline 
model  that  may  be  used  to  simulate  the  shape  of  natural  barriers  found  in 
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•    Injection  Well 
o    Production  Well 


FIGURE    19     Location  of  Combustion  Front  in  Bounded 
System 


TABLE    HI 
RECOVERY  PERFORMANCE  FOR  BOUNDED  SYSTEM 


Gas  Produced 

Oil  Produced 

Water  Produced 

Time 

(ft3) 

(bbls) 

(bbls) 

(hrs) 

0.0 

0.0 

0.0 

6.3 

0.0 

0.0 

0.0 

68.7 

19.0 

3.3 

3.3 

131.0 

113.0 

19.5 

19.7 

194.0 

223.0 

38.8 

39.5 

256.0 

334.0 

58.4 

59.7 

319.0 

453.0 

79.7 

81.6 

381.0 

556.0 

97.8 

100.0 

444.0 

668.0 

118.0 

121.0 

506.0 

767.0 

136.0 

139.0 

569.0 

858.0 

152.0 

156.0 

631.0 

930.0 

164.0 

169.0 

694.0 
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reservoirs.    By  including  boundaries,  the  area  swept  and  total  production 
may  be  determined  for  various  locations  of  the  injection  and  production 
wells. 

Comparison  with  Field  Tests 

To  establish  the  ability  of  the  streamline  model  to  determine  the 
location  of  the  combustion  front,   data  from  two  field  trials  were  used  as 
input  data  for  the  model.    The  first  test  utilized  information  provided  by 
the  Sun  Oil  Company  (17).    The  test  was  conducted  in  the  May-Libby  reser- 
voir which  has  the  characteristics  as  shown  in  Table  IV  .     The  pattern  used 
was  an  inverted  five  spot  with  dimensions  as  shown  on  Figure  22.    Of  the 
information  provided  on  Table  IV,   only  the  porosity,   injection  rate,  and  well 
pattern  and  spacing  were  utilized.    The  rationale  for  using  only  this  infor- 
mation will  be  explained  at  the  end  of  this  section.    The  saturation  and  mo- 
bility of  each  zone  was  assumed  to  be  the  same  as  that  of  the  basic  test  in 
the  first  section  of  this  chapter.    In  addition  the  amount  of  reservoir  hydro- 
carbon existing  and  consumed  in  the  propagation  of  the  combustion  front  was 
the  same  as  in  the  basic  test.    The  burned  area  is  shown  as  the  shaded 
area  in  Figure  22.      As  may  be  seen,  the  area  burned  in  the  field  test  and 
that  indicated  by  the  streamline  model  are  extremely  close. 

The  second  test  was  a  comparison  of  the  streamline  model  and  a 
field  test  conducted  by  the  Magnolia  Petroleum  Company  (9).    Table  V 
shows  the  characteristics  of  the  reservoir  and  Figure   23    shows  the  dimen- 
sions of  the  inverted  five  spot  pattern.    As  in  the  previous  test,   the  satur- 
ation and  mobility  of  each  zone  were  as  used  in  the  basic  test  as  well  as 
the  amount  of  reservoir  hydrocarbon  existing  and  consumed.    The  burned 
zones  are  quite  similar  except  in  the  region  between  the  injection  well  and 


TABLE    IV 

CHARACTERISTICS  OF  SUN  OIL 
COMPANY  TEST 


*Well  Pattern: 

Inverted  Five  Spot 

*  Well  Spacing: 

As  Shown  on  Figure  10 

*  Average  Porosity: 

31.2% 

Depth  of  Pay: 

3,400  ft 

Thickness  of  Pay: 

4.4  ft 

Average  Reservoir 
Permeability: 

1,069  md 

Average  Amount  of 
Oil  in  Place: 

740  stb/acre-ft 

*  Average  Injection 
Rate: 

8,000  cu.ft/hr 

Used  as  input  for  streamline  model 
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FIGURE    20     Streamline  Model  Burned  Zone    versus 
Sun  Oil  Company  Field  Test  (9) 


TABLE    V 


CHARACTERISTICS  OF  MAGNOLIA  PETROLEUM 
COMPANY  FIELD  TEST 


*Well  Pattern: 

Inverted  Five  Spot 

*  Well  Spacing: 

As  Shown  on  Figure  23 

*  Average  Porosity: 

27.2% 

Depth  of  Pay: 

195  ft 

Thickness  of  Pay: 

17  ft 

Average  Reservoir 
Permeability: 


Average  Amount  of 
Oil  in  Place : 


*  Average  Injection 
Rate: 


7,680  md 


1,350  bbls/acre-ft 


3,  200  scf/hr 


Used  as  input  for  streamline  model 
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FIGURE    21     Streamline  Model  Burned  Zone  versus 
Magnolia  Petroleum  Company  Field 
Test    (17) 
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production  wells    1  and   4.    The  authors  of  the  article  published  on  the  burn 
indicated  that  the  reservoir  properties  dictated  the  non-radial  formation  of 
the  burned  zone.    The  streamline  model  of  in-situ  combustion  is  based  on 
a  uniform  permeability;  therefore,   if  the  permeability  of  the  areas  in 
question  are  lower  than  the  rest  of  the  reservoir,   little  of  the  injected  air 
would  go  in  the  direction  of  wells  number    1  and   4.    In  addition,   if  the  reser- 
voir was  tighter  between  the  injection  well  and  production  wells    1    and   4 
the  production  of  these  two  wells  would  be  lower  than  that  of  the  other  wells. 
The  authors  confirmed  this  fact. 

A  third  test  was  attempted  using  the  information  published  by  the 
U.S.  Bureau  of  Mines  with  regard  to  an  in-situ  combustion  experiment  in 
oil  shale  (11,  12).    The  result  was  not  satisfactory  because  of  the  nonuniform- 
ity  of  the  permeability  of  the  oil  shale  due  to  extensive  fracturing  around  the 
injection  and  production  wells.    The  test  pattern  consisted  of  an  inverted 
five  spot  with  the  first  set  of  four  production  wells  approximately  twenty- 
three  feet  away  from  the  injection  well.    A  second  set  of  production  wells 
approximately  fifty  feet  away  from  the  injection  well  were  located  on  the 
same  axis  as  the  first  set.    With  the  extensive  fracturing  consisting  of 
electrolinking,  explosive  fracturing  and  by  hydraulic  fracturing,   areas  of 
highly  permeable  oil  shale  existed  and  enabled  excellent  communications 
between  the  injection  and  production  wells.    Since  the  permeability  in  the 
rest  of  the  reservoir  was  extremely  low,  all  injected  air  would  flow  in  the 
direction  of  the  least  resistance,   i.e. ,  toward  the  production  wells.    Making 
the  assumption  that  the  oil  shale  was  a  homogeneous,   isotropic  deposit,  as 
was  an  assumption  for  development  of  the  streamline  model,   was  therefore 
erroneous.    As  a  result,  the  test  by  the  U.S.  Bureau  of  Mines  could  not  be 
simulated  by  the  streamline  model. 
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The  articles  relating  the  two  field  trials  did  not  contain  all  the  in- 
formation necessary  to  reproduce  the  field  tests.    As  was  stated  above,  the 
porosity,  well  spacing,  and  injection  rate  were  the  only  data  taken  directly 
from  the  field  tests.    This  was  done  for  numerous  reasons.    First,  the 
streamline  model  assumed  that  the  reservoir  was  only  one  foot  thick.  Thus, 
the  combustion  front  could  be  assumed  to  be  a  vertical  front.    Although  this 
assumption  is  valid  for  a  thin  reservoir  section,  Moss  has  found  that  the 
assumption  cannot  be  made  for  a  thick  reservoir  (9).    Secondly,  the  perme- 
ability is  tied  directly  to  the  mobility  of  each  zone.    Since  this  area  is  one 
which  is  unknown  at  the  present  time,  the  permeability  of  the  reservoir  is 
of  no  use.    Lastly,  the  difference  in  production,   as  demonstrated  by  the 
comparison  of  the  Moss  and  Martin  relationships  between  rate  of  advance 
of  the  combustion  front  and  the  injected  air  flux,  would  prevent  meaningful 
comparison  since  it  is  not  known  if  the  aforementioned  relationship  is  the 
same  for  all  reservoir  materials. 

Discussion 

The  comparison  of  the  streamline  model  to  field  trials  demonstrated 
the  model,   as  it  exists,   can  approximate  the  location  of  the  combustion 
front.    Some  of  the  variables  that  affect  the  combustion  front  location  have 
been  examined  and  the  effect  determined.    The  effect  of  other  variables  that 
warrant  additional  study  are  listed  in  the  next  chapter. 

As  mentioned  earlier,  no  attempt  was  made  to  determine  the  loca- 
tion of  the  oil  bank  or  to  duplicate  production  data  due  to  the  complexity  of 
variables  necessary  to  determine  the  velocities  of  zones  2    and    3. 


CHAPTER    V 
CONCLUSIONS  AND  RECOMMENDATIONS 

Conclusions 

As  a  result  of  this  study,    the  following  conclusions  may  be 
drawn : 

1.  The  streamline  model  can  be  used  to  study  the  effects  of  vari- 
ables on  the  process  of  in-situ  combustion. 

2.  The  streamline  model  can  be  used  to  determine  the  location  of 
the  burnfront  with  an  accuracy  dependent  upon  the  accuracy  of  the  input 
data  and  the  homogenity  of  the  reservoir. 

Recommendations 

1.  The  relationship  between  the  rate  of  advance  of  the  combustion 
front  and  the  injected  air  flux  should  be  determined  for  the  specific  reser- 
voir materials  in  order  to  increase  the  accuracy  of  the  streamline  model 
for  in-situ  combustion  in  the  specific  reservoir. 

2.  A  study  of  the  problem  of  determining  the  saturation  of  various 
zones  is  necessary  to  obtain  more  realistic  values  of  the  zone  composition 
and  its  mobility. 

3.  Larger  pattern  areas  should  be  utilized  to  determine  the  best 
economic  relationship  between  injection  rate  and  pattern  area. 

4.  An  expanded  study  to  include  variation  in  reservoir  thickness 
would  be  of  great  value. 
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5.  The  effect  of  dip  should  be  considered  as  it  affects  the  flow  of 
oil  during  the  combustion  process. 

6.  The  results  of  this  study  were  based  on  a  reservoir  with  uniform 
permeability.    With  the  extensive  use  of  fracturing,   a  study  comparing  the 
effects  of  permeability  changes  is  warranted. 
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APPENDIX    A 
FLOW  CHART 


READ:   Number  of  Streamlines;    Well  Radius; 
Reservoir  Porosity,  Pressure,  Density; 
Oil  Content,    Thickness,    Boundary; 
Well  Locations  and  Rates;  Mobility, 
Saturation  of  each  Zone;    Flux  versus 
Advancement  of  Burn  Data. 


Print  Data 


Convert  English  Units 
to  Metric  Units 


t 


Compute  Pressure  at 

Injection  and  Production 

Wells 


Plot  Location  of  Wells 


Locate  Starting  Point 
of  Each  Streamline 
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Iterate  For  Each  Injection  Well 


Iterate  For  Each  Streamline 


Locate  Starting  Point 
of  Each  Zone 


Iterate  For  Each  Injection  Well 


_/\. 


1 


Iterate  For  Each  Streamline 


4     1 


4    1 
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Combustion  Front  at 
Boundary  of  Reservoir 


Yes 


Compute  Pressure  and 
Velocity  of  Standard  Fluid 


Determine  Zone  at 
Production  Well 


Compute  Conductivity  Ratio 


Compute  Air  Flux  at 
Combustion  Front 


Air  Flux  Below  Minimum 


Yes 


Compute  Velocity  of 
Combustion  Front 


4     1 


4     1 
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Compute  Velocity  of 
Zones  2,  3,  and  4 


Combustion  Front  at 
Production  Well 


Yes 


Compute  Increment  of  Movement 


Compute  New  Location 
of  Each  Zone 


Determine  Zone  at 
Production  Well 


Redefine  Production  Ratios 


Compute  Incremental  Gas, 
Oil,  and  Water  Production 


4     1 


4 

1 

2 

' 

Plot  Location  of 
Combustion  Front 

1 

Store  Zone  Locations 

i 

Iterate  For  Each  Injection  Well 

i 

Injected  Air  Produced 

< 

Compute  Time  of  Movement 

Test  Each  Streamline 

61 


Compute  and  Print  Recovery 
Performance 


APPENDIX    B 
COMPUTER  PROGRAM 


PROGRAM  COMFUInPUTiOIITPUT) 

this  program  plots  the  location  of-  thf.  burn  front 

DIMENSION  Q(7>»Xw(7),Yw<7)  »NSL (41 )  iTSUM(7 t4l)  »X  (7»M)  vY(7t4J)  , 

1  CO ( 7 1 4] ) $ C« ( 7 ? 4 1 ) • CW ( 7 1 M ) «  £y ( 7 • 41 ) • XEB ( 7 1 4 1 >  t  YEB ( 7  1 41 ) f 

2  TBl(7i41)»XFB(4l) »YFB{41) ,WATPR0(7»41) »XE2(7«41 ) «YE2(7i 

3  4U t*F2Ul) «YF2<4l) »GASPPD(7,41) »XE3<7»41) ,YE3<7,41) i 

4  XF3(41),YF3(41)iAC4(7f41)fYt.4(7»41)tXh4(4l)iYFA(<»l) 
PEAU  330i  NST»NPtNltNlMT»RI,FMR,pORiZMiPM 

READ  390»  PVOLr'X»H?XMAX,YMAX»  IPV 
READ  430,  SLOPEO»fiONEtFLOXM 
READ  440»  SLOPET»BTWO»(-L^XMN 
READ  5ro*  RHOiFUELfGALPTN 

READ  5?  6,  KATl02»RATlO3»RATlO4tWRATl2»WRATl3.*RATI4»r,RATI2tr,RATl3. 
1       GRATU 
READ  5l0«  VOL02«VOL03 

PRINT  500t  NST,NP»NI,NUATiRlfPCP,^V»PM 
PRINT  490»  PVOLMX,H,XMAXiYMAX» IPV 
PRINT  450*  SLOP£0«BOmE 
PRINT  530*  SLOptT»BTwO 
PRINT  4P0»  RHO,FUEL»GALPTN 
PRINT  470*  RAT 7 02 • RAT  103 »RAT 1 04 
READ  540»  SMBiTMH»FOPMb,FlVMB 
PRINT  3S0»  FMB,SMB»TMb»FCRMB»FlVMb 
1T0T=NP*NI 

READ  34Q»  tQ(I)iXWd)  iYW(I),  I  =  lf  i]0T) 
READ  540»  UMlN,QMINP,AREAiSDTST 
PRINT  570 

PRINT  3P0»  (Q(D  *XW(D  iYin(I)  iIsi,1T0T) 
PRINT  460»  QMIN»UMINP»A>REA»SDTST 
SVALPVsAREA»H*pOR/QMIN»3600, 
SVALUE  =  IPV°SVALPV/rJDAT 

CORRECTION  FACTORS  FROM  ENGLISH  TO  METRIC  UNITS 
GASPRD(jtK)=0.0 
WATPRD(j,K)=0.0 
MDSTX  =  MCSTX<>30.48 
MDSTYsMDSTY*30,4a 
XMAX=XMAX*30»48 
YMAX=YMAX*30.48 
SDTST"SCTSI*30,48 
h=H*30.A8 
RI=SOTST/S0.0 
NTOTaNl+fgP 
P(3s50«*30»^8 
OMlN=QMIN»7.88o55 
OMUmPsOMINP<*7,88055 
NIl=NP4l 

CTR=1.0 

PI=0.0 
PP=0,0 
00  10  1=1 t ITOT 

CORRECTION  FACTORS  FROM  ENGLISH  TO  METRIC  UNITS 
XW(I)"XW  (I  )<>30,4d 
YW(I)=YV,  (I>"30,48 
0(1)  sQ(  I')*/,  88055 
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PI=Pl-(Q(l>*ZM)/  (^^fcfc^FMB^HopM^'ALOGf  (Xw(hll).RI-X*  (I)  )«<>?♦ 
1        (YW(Ml)-YW(I)  )«»2) 

pp=PP-(Q(I  )*>ZM)/(12,5664ttFMH*»H«PM>*AL0G(  (Xw(l)*RI-xw(I))o«2* 
1       (YW(i)-Yw(I))««H) 
10   CONT  ir^)E 

CALL  HGNPLT 

CALL  PLT  (i0.0«5.5»-3> 

CALL  PLT  (-3.0,-3.00) 

CALL  PLT  (-3.0»3.0»2) 

CALL  PLT  (3,0»3.0»2) 

CALL  PLT  <3.0i-3.0»2) 

CALL  PLT  (-3.0t-3.0t2) 

SCL=. 00100^^783 

00  20  1  =  1,1  TOT 

CALL  SYMBOL  (Xw(l)«SCLfY*(l)«SCL»»21iliO.»-l) 
20    CONTINUE 

DO  40  JsNlltNToT 
NSL(J>=f\.ST*QU)/UMlN*0.45 

PRINT    370 

1B=NSL(J) 

DO    30    K=1,IR 

X{J,K)=RI<>C0S(l*6.2832oK/NSL(J)  )  **W(J) 

Y(J»K)=PI*SIN(1*6.283H*K/NSL(J)  )  *YW  (J) 

XE2(J»K)=XU»K) 

YE2(JfK)=Y  U»K) 

XE3(J»K)=X(J,K) 

YE3(J»K)=YU,K) 

XE^ (J*K) sX ( J,K) 

YEA (jtK)sY(J»K) 

OW(J,K)=o.O 

COUiKj^o.U 
CW( J,K) =0.0 

1  B  T  (  J ,  K  )  =  0  •  0 
TSUM(J,K)=0,0 

30  CONTINUE 
40  CONTINUE 
50    CONTINUE 

DO  290  j=Nll»NTOT 

1B=NSL(J) 

DO  280  KsliIB 

IF  (GW(JiK) .GT.0,0)  GO  TO  280 

XFB=X(J,K) 

YFB=Y(J,K) 

XF2=XE2(J»K> 

YF2=YE2(J»K) 

XF3=XE3<J|K) 

YF3=YE3(jfK) 

XF4=XL4(J,K) 

YF<»  =  YE4  (j,K) 

CALL  PLT  (X(J,k)°SCL,Y( J»K)*SCL»3> 

TSUMI=TSUM(J,K) 

IF  (AbS(XFB) .GT.XMAX)  270*60 
60    CONTINUE 

IF  (ABS(YFB) .GT.YMAX)  270»70 
70    CONTINUE 

VELXP=0.0 

VELYP=0.0 

VELX?=0.6 

VELY2-0.0 
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vrLX3=n.o 

VELY3=0.0 
VELXA=0.0 
VELY4=0.0 
PB=0,0 
P2=0.0 
P3=0.0 
P4=0,0 

PKSTs^M/  (l2.5664»FMBftH«Pf-») 
VKST=^M/'{6#2832*H*P0RttPM) 
00  80  L=  1  « i  TO! 

FAC4=(XF4-XW(L) >«»2* (YF4-YW(L) )*«2 
P4  =  P4-PKSTaQ(L)*AL0G(f-ACM 

VELY4  =  V[-LY<»*  ((VKST«Q(LJ/CTR)«(YF4-YW(1.))/FAC4)/P0R 
VELX^-VEL*** <<VKST*Q(L)/CTR)MXF4-XW(L> )/FACA)/POR 
VELS4SSCP1  (VELX4**2*VELY4»«2) 
FAC3s(XF3-Xw  (L)  )**2*  (YJ-3-Yw(L)  )*»2 
P3SP3-PKST*U(L)*AL0G(FAC3) 

VELX3  =  VFLX3*  <(v*ST*0(L>/CTR)«<XF3~Xw(L))/FAC3) 
VELY3=VEl_Y3*  ((VKST*Q(L}/CTR)*(YF3~YW(L>)/FAC3) 
VELS3  =  SCRT<VELx3»<>2  +  vELY3t»2> 
FAC2=tXF2-XW(L)  )  **2*  ( YF  2"Yw  (L)  >  *«2 
P2=P2-PKSTttU(L)#AL06(FAC2) 

VELX2=VELX2*  ((VKST»Q(L)/CTH)<MXF2-XV»(L))/FAC2) 
VELY2  =  VELY2*  ((vKbT<>Q(U/CTR)«UF2~YW(L>)/FAC2) 
VELS?--SfRT  (VF.LX2*»2*VLLY2*«2) 
FACbs  (XFB-Xw<L))**2*(YrB-Yw(L))»*2 
PBrPB-PKST*n(L)*ALOG(FACB) 

VELXF=VELXb*( (vKST*Q(L)/CTR)«(XFR-XW (L) )/FACB) 
VELYP  =  VELYU*  (  (VKST*Q(L)/CTR)«  (YFB-YMLM/FACH) 
80         CONTINUE 

VELSR=SQRT(VELXB*«2*VELYB**2) 

IF  (SQRf((AfjS(XF4)-PO>»«2*(ApS<YF4)-PQ)*<»2)  .LT.RI)  GO  TO  qO 
CTY={PP-Pl)/(FMB* (  ( (Pb-PD/FMR) ♦  {  (P2-PB)/Syb) ♦ (  (P3-P2)/TVB) ♦  (  (P4- 
1       P3J/F0RMB) ♦ ( (PP-P4)/FIVMp> ) ) 

GO  TO  120 
90    CONTINUE 

IF  (SORT  (  (ABS(XF3)-PQ>tt*2*  ( /IRS  <  YF  3) -PQ  )  *»2)  .LT.RIl  GO  TO  10n 
CTY=(Pp-pI )/(FMB«( ( (Pb-Pl)/FM3) ♦ ( <P2-PB)/SvB)* ( (P3-P? ) /TNB ) ♦  (  (PP- 
1       P3)/F0RMB))) 

GO  TO  120 
100   CONTINUE 

IF  (SORT  (  (ApS{XF2)-P0)«*2*  ( ARS < YF2) -PC ) *«2)  .IT. PI)  GO  TO  lip 
CTY=(PP-pI)/(FmB*(  (  (Pti-PD/FMB)  ♦  (  tP2-PB)/SMd)  ♦  (  (P3*P2)/TPB)  )  ) 

GO  TO  120 
110   CONTINUE 

DEM=FMR«( ( (PP-PB)/FMB)*( (P8-P  I ) /S^B ) ) 
CTY=(Pp-pI)/DEw 

120     continue 

velsa5vel.sb*pm/zm*cty 

IF     (VtLbA.LT»FLUXM)     140* 130 

130     continue 

velb=slopeo»velsa*bone 

GO    TO    ISO 
HO      CONTINUE 

IF    (SLOPET»vELSA.lE.ABS(BTwo )     GO    TO    270 

velb=slopet<>velsa*btwo 

lbO   CONTINUE 

AMTO 1L=G ALP TN«RHO*. 0  00  06684 
POR=POR-aMH)IL 
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FACT2  =  AMOlL<Ml.0-FUFL)*MVOLO2/  (  vO|_02*  VOL  03  )  )  /  (P0R»vnL02) 

FACT3  =  AMTOIL»(1.0-FUEL)**  (V0L03/  ( V^L02* V0L03 )  )/  (POR»vnL03)  ♦ FACT 2 

P0R=.25 

VELF^C=vEL«/VElS8 

$VELX4=VELX4 

SVELV*»  =  VELY4 

VELU  =  VELS<* 

IF  (SORT  (  (ARS(XF3)-PQ)**2*  ( ARS  I YF3) -PG ) *»2)  .LT.RI  J  GO  TO  16o 
SVELx3  =  vEL^3<'VrLFACOFACT3/P0H 
SVELY3=VELY3»VELFAC*FACT3/P0R 
VELT3  =  SCRHSVELX3«*»2*SvELY3ft*2) 
160   CONTINUE 

IF  (SQ«n(AHS(XF2)-P0)*«2*<ABS<YF2)-PQ)««2)  ,LT#PI)  GO  TO  170 
SVELx2sVEL*2*VELFAC«FACT2/P0R 
SVELY2=VELY2*iVELHAC*FACT2/P0P 
VELT2=SQRT  lSVELA2«<*2*bvELY2<><*2) 
170   CONTINUE 

5VELXBsVELX8«VELFAC/P0R 
SVELYB=VELYH*VELFAC/P0H 
VELTB=S0RTtSVELAb**2*bVELYB^o2) 
VELTN  =  VELT*» 

IF  (vtLT4.t>T.VELTM)  VELTN  =  VELU 
IF  (vELT3.<iT.VELTM)  VELTM=VELT3 
IF  (vEL.T2.^T.VELTM)  VELTr'  =  VELT2 
IF  (vELTB.bT.VFLlM)  VELTK=VELTB 

IF  (PB.^T.O.U)  GO  TO  160 
00  IPO  L  =  l»f!P 
RAD=SURT ((XFB-XW(L) )**2* (YFB-Yw (L>>«tt2) 

IF  (PAD.GT.RI)  GO  TO  1B0 
TBT (jiK)=TbUMl 
QWtJfK)rQMlN/NSL(J) 

GO  TO  2?0 
180   CONTINUE 

0T=RI/VELTM 

IF  (TSUMI*DT.GT.SVALUE)  190t20^ 
190   CONTINUE 

DT=SVALUE-!SUMI 

200   CONTINUE 

xfb=xFb«dt°svel.xb 

yfb=yfr»dt*svelyb 

XF2=XF2*nTttSVELX2 
YF2  =  YF2*DT<>SVELY2 
XF3=XF3^DTttSVELX3 
YF3=YF3*QTttSVELY3 
XF4=xF4*oTttSVELX4 

yf^hyFa+dt^svely^ 

call  plt   (*fb«scl»yfb*scl»2) 

TSUMI=TSUMl*DT 

SDTST  =  70.7*30./+0 

IP     (S<JRT  (  (A0S(XF4)-PQ><»<}2*  (ABS<YF4)-PQ)«*2)  .LT.RI)    GO    TO    210 
RATIC=0.6 

WRATlOrQ.O 

GRATlO=o,0 
GO  TO  2^0 
210   CONTINUE 

WRATI0=WRA1 14 

ratio=rati^>4 
grati0=grati4 

if   (sohm(ars(xf3)-pq)*«2*(/»bs«yf3)-pq)**2). lt.ri)   go  to  220 

GO    TO    2*0 
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220   CONTINUE 

RATlC  =  PflTlvJ3 

wRATlOrwflAl 13 

GRATI0=GRAU3 

IF  (SQR* < (ABS(XF2)-PQ)»*2* UBS  <  YF2) -PG  )  «*2) .LT.RI)  GO  TO  ?3o 
GO  TO  2^0 
230   CONTINUE 

RATI0-RATIO2 

WRATI0=WRA1 12 

GRATI0=GRATI2 
240   CONTINUE 

ADDPRO  =  CT*<C)MlN-tt(ZM/pM)/(NSL<  J)«VtLSB»POR>  >*(• 000001 9l« (VELSB*POR* 
1       CTY*>  (PM/ZM)  )♦.  O0u3*9) 

CO(J,K)=CO(J»K) ♦AD0PRO«hAtio 

GASPRD( J,K) =GA?PRD(J,K)*AOOPWO«PPATIO 

WATPRD( J,K>=WATPRD<J«K) ♦ADDPRO*WRATJO 
IF  (AHS(XFB) .GT.XMAX)  2  V  0  t  2  =  0 
250   CONTINUE 

IF  (AB5UFB)  .bT.YMAX)  2?0i260 
260   CONTINUE 

IF  (TSUHI.LTtSVALUE)  GO  TO  70 

CALL  SYMBOL  (  XFb^SCL  »  Yf- B«SCL  »  .  0  7  .  *  1 1  0  .  »-l ) 
270   CONTINUE 

X(J»K)=XFB 

Y(J.K)=YFB 

XE2(J»K)=XF2 

YE2(jiK)=Yf-2 

XE3(  JiK)rX»-3 

YE3(j»K)=Yf"3 

XE4(JtK)=X>"^ 

YE^(J»K)=Yt-4 

TSUM(J»K)*T5UMI 

CY(J.K)=CTY 
280   CONTINUE 
290   CONTINUE 

COP=C0 

CGP=0.0 

CWTRP=0.0 

CTR=0.0 

CwP=0«0 

[30  3C0  J  =  NlliNTOT 

IB=NSL(J) 

BO  300  K=1»IB 

CWP=CWP*QW(J,K)»(SVALUE^TBT (J»K) )*.  000 0062898 

COP=COP*COUtK)«.  0000062898 

CGP=CGP*GASPRD(JiK)*,0000353H 

CWTRP=CWTRP*WATPRD(J,K)*» 00000628^8 

CTR~C1R*CYU»K)/IB 
300   CONTINUE 

TIME=SvALUt/36oO.O 
CLP=COP*CwH 

PRINT  360»  CGP,C0P*CWTRP»TIME 
IF  (CWP'GT.0.0)  320t3l0 
310   CONTINUE 

SVALbL=SVALuE*SVALPV/NUAT»lo.O 
GO  TO  5U 
320   CONTINUE 
PRINT  550 
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330 
3*0 
350 


360 
370 


3b0 
390 
400 
410 
420 
430 
440 
450 

460 


470 

4b0 

490 
500 


PHI  NT 

PRINT 

PRINT 
PRINT 
PRINT 


NST,NP»NI,NUAT»RI,FI"Bi?OR»RATIOiZMiPM 

pvolmx»h,xmax»ymax, ipv 
4?6*  (Q(D  fXWn)»Vh(X)  tlsif  ITOTJ 


400* 
560 
410» 
570 


CALL  ENDPLT 

FORMAT 

FORMAT 

FORMAT 
1 
2 
3 
A 

f  ormat 

FORMAT 
1ED 


F  7.0) 

5.0tr^X,F5!0) 
CF  FLUI^  IN  ZONE 
FLUID 
FLUID 
FLUID 
FLUID 
</16X,Ll0t2ibXtLl0.2»eX»fc.l0.2il2X,E.10.2) 
(///16X»«GAS  PRODUCED       OIL  PRODUCED 
TIME  IN  H0URS«»/»16X,»cuBIC  FEET 


(4l5t4F5.2f 

( F  i  o  .  0  »  5  X  f  F 

(16X»*M0BILITY 

16X»«M0B1LITY 
l6Xi<*MOBlLlTY 
16X»*M0BILITY 
16X»*M03ILITY 


CF 
CF 
CF 
OF 


IN 
IN 

In 
in 


ZONE 
ZONE 

ZONE 
ZONE 


«Fl0.3»/» 
«Fl0.3i/i 
*Fl0.3»/» 
•Fl0.3»/i 
*F10.3) 


WATER 
RARRELS 


PRnDUC 


BARRELS**) 


FORMAT  (    /16X,FlOtli5XtF10ilt5X»F5#l) 

FORMAT  (4F7.2»I7) 

FORMAT  (   /16X»4l5,10E10.3) 

FORMAT   (/16X»5E10.3) 

FORMAT  (  /16X,e10,3i5X,E10.3.5XiE10.3> 

FORMAT  (3F10.8) 

FORMAT  (3F10.8) 

FORMAT  (//16X»17HbL0PE  LINE  0f^E  ■  F9 .7 1  3X  ,  14HLINE  CONST.  = 


Fl 


lU.B) 

FORMAT     (16X,*MINIMUM  INJECTION  RATE  =  «E10.2i/i  16v  »«MINImijm  PRO 
1DUCTION  RATE  =  *E  l  0  .  2  »  /  ♦  1  6X  »  *wFLL  PATTERN  AREA  =  «Flo.2i/»lM» 
2    *c1STaNL'E  BETWEEN  INJECTION  WELL  AND  PRODUCTION  WFLL  =  »F10,2) 

FORMAT  "(16X,*AM0UNT  OF  OIL  IN  ZCNE  2  =  «F  5  »  3  »  /  f  1  6  *  ♦  *  Af  OU\T  CF  0 
1IL  Ik  ZONE  3  =  »F5.3»/il6X»«AMCUNl  OF  OIL  IN  ZONE  4  =  »F5.3) 

FORMAT  </16X»«AVERAGE  DENSITY  UF  RES.(RHO)  =  «F  1  0  .  4  t  /  ,  1  6x  » *>  A  **0U 
INT  OF  OIL  CONSUMED  AS  FUEL  =  *F  1 0^4 »/,  1 6X  »*  AVERAGE  RES,  OIL  COMEN 
2T  GAL.  PER  TON  =  «F10«4) 

FORMAT     (/16XJORESERV0IR  PORE  VOLUME  =  *F9  ,2  »  »  /  »  \  6V  »«TP  ICkNESS  0 
IF  RESERVOIR  =  oF6.2,/tl6X,«MAXlMUN  LIST  IN  X  DIRECTION  =  «Fft.2,/, 
2   16X»#MAXIMUM  U1ST  IN  Y  DIRECTION  =  *F6.2 t / » 1 6X , « IN jECTED  PQRF  VO 
3LU^'E5  =  *>Ib) 

FORMAT  {/16X,t'NUMBER  OF  STREAM  TUBES  =  « I5»/  » 1 6X  » DUMBER  OF  PROD 
1UCTICN  WELLS  =  *l5f/tl6Xt«NUMBER  °E  INJECTION  WELLS  =  M5»/,16X»«N 
2DAT  -  «i5t/,16x»*PADlUS  CF  WELL  B^RE  IN  FEET  =  <*F5.2«/ •  16X«  oPOROSl 
3TY  OF  RESERVOIR  =  *F5#2 t/ i 16X»«coMPRESSIB ILITY  FACTO"  =  «F5,2»/il6 
4X,  ^RESERVOIR  PRESSURE  IN  ATMOSPHERES  =  «F7.0> 


F9«7«3XiUHLUE  CCNST.  =  Fl 


510 

FORMAT 

(2F10.O 

520 

FORMAT 

(9F5.3) 

530 

F  OHMAT 
10.8) 

(//16X»17HSLCPE  LINE  ]*0  - 

540 

FORMAT 

(4F]0.3) 

5b0 

FORMAT 

l/18Xf*NST     NP    NI   NCAT 

10    ZM 

PM*) 

560 

FORMAT 

l//l6X»»PV0LMX    H    XMAX 

570 

FORMAT 

<//lfiX»«  WELL  RATE 

5b0 

FORMAT 
END 

(3F1Q.4) 

RI 

YMX 
X  LOC 


FMB 
IPV 


POR 

AM8«) 
LoC  «) 


RATI 


APPENDIX    C 
NOMENCLATURE 


o 

A  Cross  sectional  area,   L 

C  Constant 

F  Factor 

h  Thickness,   L 

2 

k  Permeability,   L 

2  Length,   L 

n  Moles,  dimensionless 

P  Pressure,  M/LT2 

3 
Q  Injected  volume,   L 

o 

q  Injection  rate,   L  /T 

2      9 

R  Universal  gas  constant  (per  mole),   ML  /T  t 

SL  Slope  constant 

S  Saturation,   dimensionless 

T  Temperature,  t 

t  Time,   T 

u  Flux,   L/T 

V  Volume,   L3 

v  Velocity,  L/T 

z  Compressibility  factor,   LT  /M 

Greek  Letters 


y  Conductivity  ratio,   dimensionless 

X  Fluid  mobility,  L3T/M 
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M 

Viscosity,   M/LT 

$ 

Potential,  M/L2T2 

0 

Porosity,   dimensionless 
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Constant  =   3.  1416 

P 

Density,  M/L3 

Subscripts 

B 

Bulk 

b 

Burn 

eff 

Effective 

g 

Gas 

i 

Counter  for  number  of  zones 

k 

Counter  for  number  of  zones 

m 

Average 

n 

Maximum  number  of  zones 

res 

Value  at  reservoir  conditions 

s 

Standard  fluid 

sc 

Value  at  standard  conditions 

T 

Total 

w 

Water 

X 

x-coordinate 

y 

y-coordinate 

z 

z- coordinate 

2 

Zone  2 

3 

Zone  3 

4 

Zone  4 

5 

Zone  5 
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ADDRPO 

AMTOIL 

AREA 

BONE 

BTWO 

CO 

COP 

CTR 

CTY 

CW 

CWP 

CY 
DT 


FACB,  2,  3,  4 


Total  fluid  accumulation  for  streamline  under  inves- 
tigation,  LvT 

o 

Amount  of  oil  in  streamline  under  investigation,   L 

Surface  area  bounded  by  production  wells,   L^ 

Value  of   C    in  Equation    3.8 

Value  of   C    between    FLUXM  and  FLUXMN 

Cumulative  oil  produced  by  a  streamline,   L3 

Cumulative  oil  produced,   field-wide,   L3 

Pattern  conductivity  ratio  for  producing  well, 
dimensionless 

Conductivity  of  individual  flow  tube,  dimensionless 

Cumulative  production  of  injected  fluid 

Cumulative  production  of  injected  fluid  field-wide 
basis,   L3 

Counter  for  the  conductivity  of  individual  flow  tubes 

Time  required  to  move  one  step  on  a  streamline,   T 

Denominator  within  the  summation  of  the  velocity 
equations  for  each  zone,   L2 


FACT2,  FACT3  Term  computed  by  Equations  3.37  and  3.41, 

dimensionless 


FIVMB 
FLUXM 

FLUXMN 

FMB 

FORMB 

FUEL 

GALPTN 


Mobility  of  fluid  in  zone  5,   L3T/M 

Transition  point  of  relation  between  air  flux  and 
advancement  rate  of  burn,   L3/T 

Minimum  air  flux  required  to  move  burn,   L3/T 

Mobility  of  injected  fluid,    L3T/M 

Mobility  of  fluid  in  zone  4,  L3T/M 

Amount  of  oil  in  place  that  is  consumed  by  the  burn 

Average  gallons  of  oil  per  ton  of  reservoir  material, 
L3/M 
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GASPRD 
GRATIO 


GRATI2,  3,  4 


H 

I 

IB 

IPV 

ITOT 

J 

K 

L 

MULT 
NDAT 

NI 

Nil 

NP 

NSL 

NST 

NTOT 

PB,  2,  3,  4 

PI 

PKST 

PM 

POR 

PP 

PVOLMX 

Q 


Cumulative  gas  produced  on  a  field-wide  basis,   lVt 

Amount  of  gas  in  zone  at  production  well 

Amount  of  gas  in  each  zone 

Thickness  of  reservoir,  L 

Counter  for  number  of  wells 

Counter  for  desired  streamlines 

Number  of  pore  volumes  injected 

Number  of  image  wells 

Step  counter  along  streamline 

Streamline  counter 

Counter  for  number  of  wells 

Number  of  even  increments  of  SVALPV  at  which 
production  history  will  be  accumulated 

Number  of  injection  wells 

First  injection  well 

Number  of  production  wells 

Number  of  streamlines  emanating  from  a  well 

Number  of  streamtubes  emanating  from  a  well 

Total  number  of  real  wells 

Pressure  of  leading  edge  of  each  zone,  M/LT 

p 
Pressure  on  well  bore  of  injection  well,   M/LT 

Constant  preceding  summation  in  Equation  3.  5 

Average  reservoir  pressure,  M/LT 

Reservoir  porosity,   dimensionless 

p 

Pressure  on  well  bore  of  a  production  well,  M/LT 

Reservoir  pore  volume 
Volume  of  fluid,   L^ 


72 


QMIN 
QMINP 
QW 
RAD 

RATIO 


RATI02,  3,  4 


RHO 
RI 

SCL 
SDTST 
SLOPE O 
SLOPET 

SMB 
SVALPV 
SVALUE 
SVELXB,  2,  3,4 


SVELYB,  2,  3,  4 


TBT 

TIME 

TMB 

TSUM 

TSUMI 

VELB 

VELFAC 


Minimum  injection  rate,   L  /T 

o 

Minimum  production  rate,   L/T 

3 

Volume  of  injected  fluid,   L 

Distance  between  point  of  interest  and  production 
well,   L 

Portion  of  fluid  in  zone  at  production  well 

Portion  of  fluid  in  zones  2,  3,  4 

o 

Density  of  reservoir,   M/L 

Radius  of  well  bore,  L 

Scaling  factor  to  reduce  scale  of  well 

Distance  between  production  well  and  injection  well 

Value  of  SL  in  Equation  3.  8 

Value  of  SL  when  air  flux  is  between  FLUXM  and 
FLUXMN 

Mobility  of  fluid  in  zone  2 

(Pattern  area)*(H)*(POR) /QMIN,   T 

(IPV)*(SVALPV)/NDAT,   T 

Velocity  of  each  zone  in  x-coordinate  direction,  L/T 

Velocity  of  each  zone  in  y-coordinate  direction,   L/T 

Time  of  breakthrough  of  driving  fluid,  T 

Time  since  start  of  injection  of  recorded  production,  T 

Mobility  of  fluid  in  zone  3,   L3T/M 

Cumulative  time  since  start  of  injection,   T 

Counter  for  cumulative  time  since  injection,   T 

Velocity  of  burn,   L/T 

Ratio  of  velocity  of  burn  to  velocity  of  air  across 
the  burn 


VELSA 


Air  flux,   L/T 
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VELSAR 
VELSB,2,  3,4 

VELTB,2,  3,  4 


VELTMX 
VELXB,  2,  3,  4 


VELYB,  2,  3,4 


VKST 


V0L02,  3 

WATPRD 

WRATIO 

WRATI2,  3,  4 

X 

XEB,  2,  3,4 

XFB,  2,  3,4 

XMAX 

XW 

Y 

YEB,  2,  3,  4 

YFB,  2,  3,  4 

YMAX 

YW 

ZM 

Velocity  of  the  air  across  the  front  of  each  zone  along 
the  streamline,   L/T 

Velocity  of  the  leading  edge  of  each  front  along  the 
streamline,   L/T 

The  maximum  velocity  of  the  above  four,   L/T 

Velocity  of  each  zone  along  the  x-coordinate,   L/T 

Velocity  of  each  zone  along  the  y-coordinate,   L/T 

Constant  preceding  summation  in  Equation  3.6  and 
3.7,   LVT 

Amount  of  oil  in  zone  2  and  3 

3 

Cumulative  water  produced,   L 

Amount  of  water  in  zone  at  production  well 

Portion  of  zone  2,  3,  4  that  is  water 

X-coordinate  of  fluid  interface  along  streamline,   L 

X-coordinate  of  each  zone  along  streamline  L  that 
is  to  be  stored,   L 

X-coordinate  of  each  zone  along  streamline  L  to  be 
used  in  computations,  L 

X-coordinate  of  maximum  point  of  interest,   L 

X-coordinate  of  well  locations,   L 

Y-coordinate  of  streamline  L  on  well  bore,   L 

Y-coordinate  of  each  zone  along  streamline  L  to 
be  stored,   L 

Y-coordinate  of  each  zone  along  streamline  L  to 
be  used  in  computation,   L 

Y-coordinate  of  maximum  point  of  interest,   L 

Y-coordinate  of  well  locations,   L 

Compressibility  factor  of  driving  fluid,   LT  /M 
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